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This paper gives a constructive answer to the question whether photon states can 
. contain or not, and to what extent, the readings of rulers and clocks. The paper first 
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shows explicitly that, along with the momentum representation, there is room in the 
one photon Hilbert space for an alternative position representation. This is made 
possible by the existence of a self-adjoint, involutive, position operator conjugate 
to the momentum operator [1]. Position and momenta are shown to satisfy the 
Heisenberg-Weyl quantization rules in the helicity basis, which is analyzed anew 
from this point of view. The paper then turns to the photon's time of arrival. By 
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, picking an appropriate photon Hamiltonian - using Maxwell equations as the photon 



Schrodinger equation - a conjugate time of arrival operator is built. Its interpretation, 
including the probability densities for the instant of arrival (at arbitrary points of 
3-D space) of photon states with different helicities coming from arbitrary places, is 
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Introduction 

Since the advent of Special Relativity, light propagation has played a very special role 
among the variety of natural phenomena studied by Physics. In conjunction with rods 
and clocks, light signals were a main component in the construction of the -then- new 
theory. The universality of c, the speed of light, set the new standard commensurate to the 
breadth of relativistic physics. While light was understood at the time as a manifestation of 
electromagnetic waves, the interaction of matter and radiation demanded a quantum leap: 
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the photon. This is the place to recall the mismatch between the prominent role played 
in relativity by the arrival times of light signals and the position of their group centers 
and wave fronts, and the lack of the corresponding operators in quantum physics [2, 3, 4]. 
The reason of this deficiency of the quantum theory [5] can be traced back to the fact that 
the photon does not have a complete set of helicity states; the state with helicity A = is 
missing. This paper aims to give a constructive answer to the long debated question whether 
the photon state can contain or no, and to what extent, the readings of rulers and clocks, 
i.e. the where's and when's of the photon. In view of the recent experimental advances, an 
appropriate formalism is necessary for a meaningful discussion of the propagation and arrival 
of light signals, mainly if these arrive photon by photon. The literature on this question is 
abundant and only a part of it is in the bibliography at the end of this paper. The reader 
is referred to [6] for a concise and clear introduction to this topic. 

Electromagnetic waves propagate at the speed of light in vacuum. This basic tenet 
of relativity, has been challenged recently both from the experimental and the theoretical 
front. It is remarkable that the theoretical part of this question was settled long time ago 
by Sommerfeld [7] and Brillouin [8] , who demonstrated that the electromagnetic wave fronts 
always move with speed c. To my knowledge, the cases where the group velocity in vacuum 
was found larger than c where always produced by an incorrect identification of the signal 
center. Superluminal propagation in vacuum turned always to be a geometric artifact of the 
wave form and not the group speed of the wave packet. The main point to highlight here 
is the absence of foundational problems associated to the definitions of the group position 
and time of arrival of classical e.m. waves. It may be difficult to ascertain these questions 
in a given situation, but only from the technical point of view, not as a matter of principle. 

Electromagnetic waves are composed of elementary quanta: photons. Like the rest of 
elementary particles, these are mathematically identified as elementary systems under the 
Poincare group [9, 10]. Namely, all the possible sets of values that the dynamical properties 
of the system can attain are connected by transformations of the group to an arbitrary, fixed, 
set of values, taken as standard. The different classes of elementary systems are identified 
by the values taken by the Casimir operators of the group. Photons are characterized by 
lightlike momenta p 2 = (p ) 2 — p 2 = and a helicity value that is either A = +lorA = — 1, 
there are thus two irreducible representations (1,0) and (0,1) (in the standard notation 
(j,f)), that are combined into an irreducible entity by parity (1,0)— > (0,1). Finally, photons 
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transform under the representation (1, 0) + (0, 1) of the complete Poincare group. 

The properties of photons are thus given by the eigenvalues of the subsets of commuting 
generators of the Poincare group in that representation. There are two widely used alterna- 
tive subsets: a), {p, A}, that describe plane waves, suited for the description of initial and 
final scattering states, free flying photons etc, and b). {p°,j,j z , parity}, that describe spher- 
ical photon waves that simplify the analysis of radiation, and in general, of those systems 
where there is a singled out space point. 

The rapid pace and depth of the current advances in the manipulation of photon states 
are demanding a parallel progress on the theoretical side. In some cases it is necessary to 
describe accurately the behavior of photons that are in the near region. In others, their 
space time distributions in regions close to the microscopic realm, for states containing just 
one photon or a few of them, in states possibly entangled that spread over these regions, etc. 
Current experimental set ups include a variety of very, sensitive single photon devices. Even 
more, many of the most exciting results have been obtained experimenting with these highly 
nonclassical light states. On its side, the theoretical framework should be able to address 
simple questions like the time of arrival of a photon at the place where some device is located, 
the relation between the orbital angular momentum of the photon and its position at some 
time, etc. It is necessary at least to give a meaning to these questions, something whose 
very possibility has been the subject of debate [2, 11, 12, 13, 14, 15] since the early times 
of Quantum Mechanics. I shall avoid here repeating the well known pros and cons for the 
wave function of the photon [16, 17, 18, 19, 20, 21] and whether it represents the probability 
amplitude for finding the photon in a given region of space (technically a Borel set A e -R 3 ). 
Instead, I shall trace the problem back by an alternative route to non-relativistic quantum 
mechanics (NRQM) where the vanishing photon mass (M = 0) prevents the definition of a 
position operator q for the photon. 

Let me first give some support to the above assertion. The dynamical properties of an 
elementary object, i.e. a free particle, in non-relativistic quantum theory can be traced 
back to its behavior under the Galileo group [22]. In particular the position operator q 
is related to the Galileo boost operator G by Gi = Mcji, where M is the mass of the 
particle. For sure, it can be defined directly by qi\x) = Xi\x), so that its components are 
commuting, its eigenvalues real numbers, and its eigenfunctions Dirac delta functions. On 
the other hand, this definition has to be compatible with the fact that the momentum 
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operator p generates translations in configuration space, [qi,pj] = ih. Accordingly in the 
momentum representation, more appropriate for a discussion of the Galileo group generators, 
the position operator acts as a derivative on the momenta = ikd/dpf, if pi is the basic, 
multiplicative operator of the representation, then fa is derived from it. In other words, 
particle position emerges from the Galilean properties of the particle, specifically from its 
behavior under space translations. 

As a consequence of the above, one may wonder that if in NRQM a massless system can 
not have a properly defined position operator, then: Why expect the opposite in RQM? The 
contenders in the debate about the photon wave function could be aware or not about the 
lack of a non-relativistic q for the photon. In any case, they were undoubtedly influenced 
by their consequences, and induced to give up in the search of a relativistic q for M — 0. 
The fact that the photon little group is E{2) instead of SU{2) added the filling drop to the 
confusion. In some studies of the position operator for relativistic systems, photon transver- 
sality was identified as the obstacle to the existence of a photon position operator [4, 5]. 
Other analysis [3, 20] pointed to the lack of commutativity of its components (or better, of 
the Lorentz boost components) as the reason for the lack of an appropriate q. Finally, M. 
Hawton found out [1, 23, 24] very recently a self-adjoint position operator with the required 
commutation rules, that operates in the photon Hilbert space as a physically acceptable 
position should do. This opens up the possibility of analyzing the question of photon lo- 
calizability from a new and more powerful setting: A good position operator should come 
endowed with a probabilistic interpretation, which is the tool suited to the analysis of lo- 
calizability [5, 13, 14, 15]. On the other hand, the most recent theoretical analysis [25, 26] 
indicate that individual photons can be localized in regions much smaller than previously 
thought, with entanglement (instead of wavelength) as the physical ruler [27] in the exper- 
iment [28]. This increases the value of the Hawton construct and the urge to explore its 
consequences. 

In this paper I shall analyze the properties of the photon states in Minkowski space-time. 
On the way, we shall devise several tools necessary to this end. The photon Hilbert space will 
be introduced in Sect. II, where we also present two alternative conjugate representations 
and show the probabilistic interpretation that they can be given. Sect. Ill is devoted 
to the explicit construction of the Hawton position operator in these representations. We 
revise the properties of the operator in the helicity basis and conclude that it qualifies as 
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a good position operator. Finally, we complete the standard picture of one photon states 
with additional position dependent information. Conversely, Sect. IV is devoted to the 
time dependent information. We first consider Maxwell equations as the photon Schodinger 
equation and solve the time evolution associated to it. Then, we set up the formalism for 
analyzing the photon time of arrival at an arbitrary point of space. It will be necessary 
to split the Hilbert space into two subspaces putting the eventually detected states in one 
of them. In it we build the sought time operator, obtain its eigenfunctions and give the 
positive operator valued measure that permits a probabilistic analysis of the times of arrival 
of a photon in a given state at an arbitrary space position. The paper ends in Sect. V with 
a summary of the results, some considerations on their meaning, and indications for the use 
of the formalism in the analysis of some current experiments. 



One photon Hilbert space 

We shall work in the Coulomb gauge i.e. within the Hilbert space H of one particle 
transverse states A^p), i = 1,2,3 defined on the forward light cone (p° = \p\) (note the 
symbol", it is a label for transversality) . The scalar product in Ti is defined as 

(A,A')= J da(p)A*(p)A' t (p) (1) 

where da(p) = d 3 p/2\p\ is the measure on the light cone. An arbitrary photon state A E H 
can be written as 

Mx) = (2tt)- 3 / 2 / -^L e *» Mp) (2) 
J ^p\ 

Note that A^x) is solenoidal so that, in spite of being (A, A') = J d 3 xA*(x)A'i(x), x does 
not qualify for a good position (in poor words: XiA(x) is not solenoidal, so x« takes A out 
of the Hilbert space). 

Due to Poincare invariance (that we are taking fully into account, in spite of the non 
covariant notation), the helicity W is a good quantum number. In fact, p and W = S(p/\p\) 
- where Sfj = —ie ai j are the spin-1 matrices - form a complete set of commuting operators 
on the mass shell. In the following, I shall work in momentum space where the momentum 
operator p is simply represented by the vector variable p = (pi, p 2 , Pz) ■ An arbitrary operator 
G function of the momentum, that commutes with the helicity [G, W\ = 0, shall have 
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eigenfunctions V GX (p) satisfying 

(P a v G ,,y(p) = p a v^ x (p), (wv G , x y(p) = xv^(p), (G%) i (p) = cr G » (3) 

The second of these equations implies that V GX (p) oc e l (p, A), where e(p, +1) (e(p, — 1)) 
are the right (left) polarization vectors for momentum p. There is some arbitrariness in 
their definition, the convention of Ref. [29] adopted here is e;(p, A) = Djj[A p< _ fe ] e,-(A), for 
p = Ak and = (0, 0, |A _1 j>|), with e(A) such that £3 e(A) = Ae(A). In spherical coordinates 
where p = (k, 9, <p) and k = (0, 0, k) - i.e. taking A as a pure rotation for simplicity - the 
polarization vectors are given as linear combinations of the unitary vectors e(p, a), a = 
9, ip, k (or cr = l, 2, 3): 

e(p,k) = (dp/dk), e(p,9) = (1/jfc) (dp/09), e(p,<p) = (l/ksin9) (dp/dip) 

in the standard form: 

e i (p,\) = --±={e i (p,9)+i\e i (p,<p)} A = ±l (4) 

The representations of the Poincare group for the photon states are the (1,0) and (0,1), 
that correspond to A = — 1 (left polarization) and A = 1 (right polarization) respectively. 
They are combined in the direct sum (1, 0) © (0, 1) that becomes irreducible when parity is 
included in the group. No other states exist (see for instance [29] page 69) because the little 
group is not semi-simple. In particular, the polarizations can not be taken as four-vectors. 
Even if we form the object e M (p, A) = (e°(p, A), e(p, A)), it not only has to satisfy in every 
reference frame p^ l e fl (p,\) = 0, but also e°(p, A) = 0, or any other condition chosen to fix 
the gauge. In short: in spite of its four-dimensional appearance, the polarizations do not 
transform as four- vectors under Lorentz transformations, but as a sort of connections [30, 31]: 

A£e"(Ap, A) = exp{iA©(p, A)}e^(p, A) A; A) (5) 

where Q(p, A) is the Wigner rotation angle, and il(p, A; A) is a gauge transformation fixed 
by the gauge condition (e° = 0). Of course, the photons can be given a covariant (ten- 
sor) representation by means of the gauge invariant antisymmetric tensor F a p(p, A) oc 
i^ a pnu(p^^(p,^) — P u ^(p, A)), that gives rise to the electric and magnetic fields. With 
these caveats in mind, I shall continue to use the three dimensional notation to work in the 
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Coulomb gauge in the following, turning to [30, 31] for questions of Poincare covariance or 
gauge invariance. 
By construction, 

(S • e(p, k))ij = (S k ) i:j , and (S k )ij€j(p, A) = Ae;(p, A), A = ±1 (6) 

We can now write Vq X (p) = gcip) £ l (p, A), where gcip) is a G dependent function of p to be 
determined in each case. By these definitions, these functions span the transverse subspace 
orthogonal to p, and are eigenf unctions of the helicity with eigenvalue A. 

Momentum and position representations 

In the trivial case of the momentum operator G = p where Vp X (p') = g p {p') e*(p', A), it 
is straightforward to obtain g p . We want (Vp^, Vp 2i A 2 ) = ^MM^ip^ ~~ P2), but we only 
have 

iY PlM ,V P2M ) = J da(p)V;i Xl (p)V; 2M (p) = 5 XlX2 J da(p)g; i (p)g P2 (p) (7) 

where we used that e*(p, A)e(p, A') = 5\\r. This requires the value g p (p') = \^2\p'\5^ (p'—p) 
for g p , so that, apart from phases, the eigenf unctions have to be 

VUp ) = fiW^P', A) 5 (3) (P' - p), and V^x) = e\p, A) e *» (8) 

Note that the second equation in (8) comes from the first one by straight application of Eq 
(2). 

The next step is to assume the existence of a position, namely a vector operator G = 
q = (<fi, ?2, 53), with simultaneous eigenfunctions in 7i that can be given as V^ x {p) = 
9q(p) A). For this to be possible it is necessary that: 

i. Position and helicity commute, that is [q\ W] = 0. 

ii. The three components are in involution, that is [<f , = 0. 

To these conditions we have to add that positions and momenta be canonically conjugate 
operators: [q 1 ,^] = iS %3 \ Finally, it will be necessary to check whether q gives rise to a 
probabilistic interpretation or not. 

We shall return to the explicit construction of the position operator later on but, for 
the time being, we assume its eigenfunctions exist in the transverse Hilbert space of the 
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photon and form an orthogonal and complete set. We first explore the implications of the 
orthogonality. A direct computation using V* x {p) = g q (p) e' l (p, \) gives 

(VquM, V Q2i x 2 ) = 5 Xl x 2 J da(p)g* Qi (p)g q2 (p) (9) 

but, as in the case of the momentum operator, we would like to have (V qit \ 1 ,V q2! \ 2 ) = 
8\ 1 \ 2 $ ( ' 3 Hqi — Q2)- This requires g q {p) = (27r)~ 3 / 2 \j2\p\ e~ tpq . We can now give the position 
eigenfunctions in both, momentum and coordinate representations: 

V* A (p) = (2vr)- 3 / 2 V^jif e*(p, A) e-** and V^x) = (27r)~ 3 J d 3 pe*(p, A) e*-* (10) 

It is clear from the last equation that x does not correspond to the position eigenvalue q. 
This mismatch is due to the presence in the integrand of the rhs of momentum dependent 
polarization vectors that link momentum with spin. This equation summarizes much of the 
troublesome nature of the coordinate space wave function of the photon. 

The above eigenfunctions form complete sets. By straightforward calculation using 
J2\ e"(p, A)e J (p, A) = 5 tj — p' l p> /\p\ 2 = 5 l l(p), we get the decompositions of the identity 

ncpi.p*) = £a/^(p)^>i)^> 2 ) = ^'(Pi)5 (3) bi -P2) (ii) 

A(Pi»P2)= ZxS&qVitipJV&fa) =2\p l \5 l l( Pl )5^(p l -p 2 ) (12) 

We have arrived at two alternative complete sets of commuting operators (momentum, 
helicity) and (position, helicity), whose simultaneous eigenfunctions form alternative bases 
of the Hilbert space Ti of the massless (1,0) © (0, 1) representation of the Poincare group. 
We recall again the difference between position and coordinate, and refer the reader back 
to Eq.(10) for a check. Given an arbitrary state \& in Ti, we could employ Dirac bracket 
notation to denote its components in both representations 

(pA|¥) = J da(p')V;* x (p')¥(p') = _^*( p ,A)^(p) (13) 

(q\\V) = J d<j(p')V£ x (p>)¥(pi) = (2vr)- 3 / 2 | d 3 p (pA|¥) (14) 

The relation (14) is the standard one between position and momentum representations in 
quantum mechanics. In fact, 

(qX\p\') = (2tt)- 3 / 2 5\\r exp(ipg), (pX\p'X') = 6^6® (p - p'), (q\\q'X') = 5 xx ,5^(q - q') 

(15) 
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Note also the scalar product in both representations: 

<V#> = E / d 3 p\( P X\n 2 = E / d 3 q\(q\\n 2 (16) 

A A 

The absence of polarization vectors in the integrand in the rhs of (14) eliminates the ob- 
structions [5, 13, 14] that prevent from considering 

P*(A,A)= / d 3 q\(qm\ 2 (17) 

J A 

as the probability for finding a photon of helicity A in the Borel set A e R 3 . In particular, 
P</,(A) = is possible. The reason is that, in the absence of polarization vectors in (14), 
the Fourier transform of (q\\^/) can be an entire function of p. Then, according to the 
Plancherel-Polya theorem [32], P^,(A) may vanish in Borel sets like A e R 3 . Physically this 
is necessary in order to cope with the absence of photons in A. It is the product of ^ l (p) 
and e l (p, \)/^2\p\ in (13) that has to be entire. This can be so even if the presence of 
\J~\p\ explicitly and in the e's implies that the individual &(p) can not be entire functions 
and, therefore, can not be given a probabilistic interpretation. In any case, it is worth to 
recall that the expressions in (13) are for fixed helicity. Thence photon states with fixed 
helicity could be localizable in spite of [13] and of theorems 1 and 2 in [14]. This is a first 
consequence of the until now only hypothesized existence of q. 

The expansion of arbitrary states in Ti in terms of momentum and helicity eigenstates 
can be given inverting (13) (recall that H is the space of transverse states) 

* i (p) = V^|plE c< (p. A )<pW ( 18 ) 

A 

a relation that is customarily used in conjunction with Eq. (2) to write 

**(a!) = (2tt)- 3/2 E / d 3 p e ipa V(p,A) (pA|tf) (19) 

A J 

This is the standard coordinate representation used to describe the one photon states and - 
promoting (pA|\I/) to the ranks of annihilation operators - the photon field. The comparison 
of (14) and (19) indicates clearly that it is the presence of the p-dependent polarization 
vector in the integrands of (10) and (19) that forestalls the interpretation of a? as a true 
position. The decompositions Eqs. (18) and (19) clearly show that the Sb(p) are transverse 
(i.e. J2iPi^i(p) — 0) and the ^(a?) solenoidal (i.e. d^> i(x) / dxi) = 0. The interested 
reader is referred to [20] for the use of (19) as a photon wave function and to Landau and 
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Peierls [12] to explore the meaning of the non local relations between 4f(x) and (q\\^). It 
is straightforward to show that they are 

ft(x) = E/ d 3 qV* x (x)(qm, and (q\\^} = J <PxV£(x)ft(x) (20) 



Position operator 

In the previous section we have seen that there is room in the Hilbert space of one photon 
states for an ordinary position operator. In fact 

{q , m\^) = q , l {q , m, and (pA|^> = ih(-^)(p\\rl>) (21) 

Here we will analyze the features of this operator as seen in the space of transverse states. 
Following Hawton [1], we first identify its structure by letting it to operate on the eigen- 
functions V q \. Then we shall show that it qualifies as an appropriate position operator. 
Concisely, our first task is to find the operator q a , a = 1,2,3 whose eigenf unctions are the 
V's given in Eq. (10) 

{q a V q , x } i (p) = q a {V,, A }.(p) (22) 
By using above the explicit form (10) of V Qt \, Hawton got after some algebra [1] 

{q a V q ,x} i (p) = V + (Q<%) {V q , x } . (p) =► {q%j = i5 tJ V a + (Q«%- (23) 

where 

V a = ^2\p~\^^=, and(Q a ) l3 =t ]T e,(p, a) {V%(p, a)} (24) 
°P \ 2\p\ a=e,<p,k 



As is well known, the operator ordering implied by V a is necessary to make it self-adjoint 
with the measure da(p). Notice that the attained operator (23) is independent of the helicity 
quantum number, being a matrix in the coordinate indices. By computing the derivatives 
of the basis vectors that appear in the definition of Q a , one obtains the explicit expression 
of the operator 

(q a h = iSijV a + -L [e(p, k) A S l3 ] a - ^-e a (p, V )W l3 (25) 

IpI IpI 

due to Hawton [1, 23, 24], who correctly identified the first two terms as the Pryce posi- 
tion operator [3], and the last one as a compensating term [26] for the topological photon 
phase [33, 34]. She also realized that, due to the last term on the rhs of (25), q appeared as 
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a set of three components in involution, something that previous position operators did not 
meet. 

A compact expression that summarizes much of the above, shows explicitly the relation 
between this position operator and the spinless one iV a , and gives a rationale for it, is 

(q a h=* E ei(P,<r)V a e,(P,<r) (26) 

<T=0,(fi,k 

The sum in (26) includes the longitudinal polarization e(p, k) but, in spite of it, q operates 
within the transverse subspace Ti overcoming some queries put forward by Wightman [5]. 
Using (26) on an arbitrary function \I> G H, expanded according to (18), we obtain 

{<f*}>) = E/ d 3 P'KAP)^(P'^\^) = >AWE^A)<^(p,A|*> (27) 

which shows explicitly the transversality of q ■ We also recall here that [q a , W\ = as 
can be seen by explicit computation using (23) and the definition of the helicity. Putting 
all together, it is possible to use the familiar notations: 

<f = £ [d 3 p\p\)i—(p : \l and q a = J2 [ d 3 q \q\)q a (q,\\ (28) 

A °P X 

for the position operator in the helicity representation q a . Finally, the wave functions can be 
interpreted as in the non relativistic case: given a photon in the state Py(p, A) = | (p\\^) | 2 
gives the probability of finding it with helicity A and momentum p, and Py(q, A) = | (q\\^f) | 2 
gives the probability of finding it with helicity A in the position q. 

Any arbitrary operator O defined in the Hilbert space H can be represented in the two 
different bases introduced above. Using the scalar product (1) and the relation (18) we get: 

*> = (2tt)- 3 / 2 / da{p) J2 ®*(P)6^ 3 (P) (29) 

ij 

= (2vr)- 3 / 2 / rf 3 pE(*IP A ) Oxy{pXm (30) 

so that there is a well defined relation between the operator's expressions in both bases: 

d xy (p,d/dp) = ^=^ e *(p,A)a i ,(p,9/9p) ej (p,A') V / 2W (31) 
V 2 IpI a 

Applying this to the canonically conjugate momenta = 5ijp a and position (26) operators, 
we get them in the helicity basis as 

P a xx> = hx'P a ; and q a xx> = (32) 
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We have thus recovered the old Heisenberg-Weyl quantization rules as anticipated in (21). 
They are valid in the helicity basis and only in it because, due to helicity conservation, addi- 
tional terms appear in other frames to compensate for the effects of d/dp on the momentum 
dependent polarizations. As said above, the meaning of the derivative in the helicity basis 
is that of a covariant derivative. In fact, it is related to the standard covariant derivative 
introduced in ref [26] to account for the topological phase [33, 34] of the photon. 

d 



XX' ~ °AA' 



dp a 



+ XD a (p) 



(33) 



where D a (p) = cot 9e a (p, <p)/\p\. Notice, by the way, that the helicity operator Wij trans- 
forms to the helicity quantum number A in the transformation from the spin basis onto the 
helicity basis. The spin matrix also undergoes a similar transformation 

S- xy = 5 X y\e a (p,k) (34) 

The transformation to the helicity basis analyzed above — > 0\x' should not be mistaken 
with the similarity transformation to the photon frame, namely: 

6 a aa , = £ ei (p, <7)d?- ej (p, a'), with {a, a'} = {k, 9, tp} (35) 

ij 

This is simply the rotation from the fixed cartesian axes to the axes lying along the unit 
vectors e(p,a). A most striking case occurs for the spin matrix that, by (35), becomes 
S^a' =i T,a" <w a" e a (p, o") ■ Notice the difference with (34): 

'S a , S b ] , = i £ e abc S c aa , while \S a , S b ]^ r = (36) 

J (7(7 L J A A 

c 

The vanishing of the second commutator is of no surprise as there are no remains of the 
spin matrices in the helicity representation. The dimension of the spin space is 3, while the 
helicity basis is a sum of two one-dimensional representations. Even after putting together 
the two parity related representations (1,0) and (0,1) to have both helicities, the would- 
be helicity eigenstate is outside the representation space. The relation of these facts 
with gauge invariance was discussed [30, 31] a long time ago. From the geometric point of 
view, the transformation Oij — > Ow is a projection from R 3 onto the transverse subspace 
spanned by e(p,9) and e(p,<p), along with the appropriate label rearrangements. As a 
result, the operator products must be handled with caution: in general their transform shall 
not coincide with the product of the operators' transforms. We saw this in Eq. (36) for the 
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spin matrices. The same occurs for the angular momentum and the boost operators and, in 
general, for all the operators that pull the states out of the transverse space. This sometimes 
led to define the helicity representation through a non singular transformation applied to 
the spin representation: 

<%(p, d/dp) = — L= £ R- 1 (p)d rs (p, d/dp)R v (p)^2\p~\ (37) 

J2\p\ rs 

where R(p) can be any of the rotations from the standard momentum (0,0, |p|) to p. De- 
noting by e(a), a — 1, 2, 3 three unitary vectors along the fixed axis 

ei(p,a) = R ij (p)e j (a) (38) 

Notice that the intrinsic arbitrariness in R, due to the invariance of the standard momentum 
under rotations around it, is removed by the choice of fixed vectors ej(cr) in (38). Note also 
that the helicity representation operators are obtained by projecting (37) on the fixed helicity 
basis: 

6xx> = E <W 0' i3 6,(A') (39) 

ij 

where, according to (4) and (38), e;(A) = -^{ej(l) +iAej(2)}. Needless to say that (39) is 
invertible only within the subspace orthogonal to the standard momentum. 

The application of the above results to the specific case of the electromagnetic field 
completes the standard picture of the one photon state with additional, position dependent, 
information. A photon in a state A can be given as 

A(x) = (2vr)- 3 / 2 £ / -^E= e*» e\p, A) (p\\ A ) (40) 
or in the alternative form 

A(x)=Y, fdqV^x) (q\\ A) (41) 

A J 

this reinforces the interpretation of q as a position operator and of V q \{x) (given explicitly 
in (10) as the configuration space amplitude of the photon state localized at q with helicity 
A. 

Note that we could add zero components to the polarization vectors to form a fourvector- 
like object A^(x). However, as shown in (5) this object does not transform as a fourvector, 
but inhomogeneously as a connection: 

U[A]Ap(x, XMA- 1 } = K (2n)-^J e^ Ax {e u (p,X) - Vv Q(p, A; A)} (pA| A) (42) 
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A tensor-like covariant object for the photon can be defined as 

F^x) = (27r)- 3 / 2 ]T / -p= e*»i {p, e„(p, A) - p v e,(p, A)} (p\\ A) (43) 

A J yj2\p\ 

This contains the same information that A^, but transforms as an antisymmetric tensor; its 
components constitute the electric and magnetic fields. In the next section we shall use (43) 
to construct and solve the photon Maxwell equations in vacuum. 

Time operator 

The Galileo boost operator in the NR representations for elementary systems has only 
three components Gi, i = 1,2, 3. Each of them is associated to a component of the position 
operator. There is no room in the representation space for an additional boost associated 
to the time. In fact, time is invariant under Galileo boosts. Most likely, this is hidden 
among the reasons [2, 35, 36, 37, 38, 39, 40, 41, 42, 43] for the difficulty of finding a time 
operator that behave just like the other position operators, something tempting in the study 
of individual particles. As explained in [43] , this is a misleading approach to the role of time 
in quantum mechanics and its use led to several dead ends in the past. A question to be 
taken into account is that time translations are given in terms the Hamiltonian H which 
is not independent, but is fixed as a function of the other operators in the representation 
by the mass shell condition. The same is to be expected for a time operator, instead of 
the addition of a new independent component to the position, a given function of the basic 
operators, whose explicit form depends on the very properties of the representation. 

Leaving aside the far-reaching but still out of reach question of the status of Newtonian 
time in quantum mechanics, and the search of its corresponding operator - if any - it is 
possible to identify [44, 45] time-like properties of quantum systems. The simplest of these 
appears for free elementary particles in the form of the time of arrival at a fixed position 
X, a deceptively simple property whose analysis in quantum mechanics is very delicate 
and full of traps. Classically it is a derived quantity, a function of the initial position and 
momentum whose values are either "never" (when X is out of the particle trajectory), or a 
real number that solves the equations of motion for t as a function of X. Notice that some 
kind of integrability is implicit for this classical notion to work properly [46]. Quantization 
requires the splitting of the Hilbert space into two orthogonal subspaces, one that contains 
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the never detected states, and other that contains the eventually detected ones. The time of 
arrival acts within this last subspace [44, 47], where it is represented by an operator which 
is not self-adjoint but only maximally symmetric. Its eigenstates are not orthogonal, so that 
instead a projector valued measure, the statistics of the times of arrival at a detector can 
only be analyzed in terms of positive operators [48] . Below, I shall show how this works for 
the photon. 

First, we shall study the time dependence hidden in the bracket (pA| A(t)). Taking into 
account that p 2 = 0, and that pe(p,X) = 0, we get from (43) the equations d^F^ = 0. 
These cover the full set of Maxwell equations due to the duality relation among null fields 
in vacuum [31]: 



Fij — tijk{Ek + iX B k ) 

(44) 

F ok = -iX(E k + iXB k ) 

We now have two equations d l F i0 = and d°F 0j + d^F^ = 0. Calling F(x, X) = E(x, X) + 
iXB(x, A), the first equation reduces to the divergenceless condition 'VF(x, A) = and the 
second to the Schrodinger equation [3, 16, 18, 20, 21]: 

i dF ^ X) = AcV A F(x, A) (45) 

Recalling the expression for the spin matrix (S^)i k = i€ij k , this equation can be given a more 
transparent form: 

= HSphFj(xA) (46) 

where p = —iV x . Using now (43) in the definition of F we obtain 

F(x,X) = (2tt)- 3 /2| d 3 p^2\p~\ze^e(p,X) (pX\A(t)) (47) 

therefore, (Sp)ijFj = A|p|Fj and hence, 

l ^^ = c\p\F(x,X) (48) 

Therefore, the photon Hamiltonian is H = c\p\5\x' in the helicity basis. 

Due to the simple expression of the photon Hamiltonian, the position operator in the in 
the helicity basis evolves quite simply: 

= i\q\cJp^} = c ^, ^! = q a (t) = q a + c f^t, p a (t) = p a (49) 

at v \p\ at \p\ 
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In the first of these Eqs. we displayed explicitly the hamiltonian obtained from the mass 
shell condition H(p) = c^fp 1 , but then we continue to denote it by \p\ as before. q a and p a 
are the position and momentum operators of the photon at t — 0, while q a (t) is the position 
operator at a time t (we are working in the Heisenberg picture). All this is trivial and 
prompts to the seemingly innocent question: When does the photon arrives at a position xl 
This is the alternative to the standard quantum mechanical question: Where is the photon 
at time t? Much in the same way that the question "where?" brings a position operator 
whose probability distribution is used to ascertain at what place!, the question "when?" 
demands the introduction of a time operator, with the meaning of time of arrival (ie. at 
what instant 1 ?), to close the logical loop. Several obstructions prevent the existence of such 
an operator in the standard formulation of quantum mechanics. The earliest one, which 
turned out to be the strongest, was formulated in the early days of quantum theory by 
Pauli [2]: The simultaneous existence of unitary representations for both, energy and time 
conjugate operators, is incompatible with a bounded hamiltonian, opening up at least a sort 
of infrared catastrophe (see however [43]). This led to renounce to the self-adjointness of 
the time operator, keeping its maximal symmetry only (i.e. T* T —T). In these conditions, 
time eigenstates are no longer orthogonal (unless some steering regularization is used [44]) 
and, instead of the traditional projector valued measure, it is necessary to turn to a positive 
operator valued measure [48] to interpret the formalism. The interested reader is referred to 
the review [49] to get an account of these and related issues, and references to the original 
literature. 

Generally, the search of a time of arrival operator has been undertaken in the case of 
one space dimension. The present author studied the case of the free particle in three space 
dimensions [47] concluding that the detected states are confined within a subspace of the 
whole Hilbert space. Outside it is the realm of no detection, that is, of those states for which 
the time of arrival is "never." In classical terms, they miss the detector, whose efficiency 
- a different question - is assumed to be 1. On physical grounds, detection requires a 
constraint: that the particle momentum is parallel to the line joining q with the arrival 
(detector) position z. This vector constraint and the free particle Hamiltonian form a first 
class system. The use of this formulation made possible the obtention in [47] of a very simple 
solution for the time of arrival in three dimensional space, basically an extension of the 1-D 
results. 
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A time operator for photons requires of a 3-D setting. Not only because the transverse- 
vector character of the electromagnetic field stripes the 1-D approach of credibility. Also from 
the practical side: we shall analyze later on the arrival of photons through inhomogeneous 
media, where direction changes will in general take place, whose description needs of more 
than the mere distance covered. Therefore, we shall examine the first class constrained 
system that evolves in the detected subspace as the first step towards the time of arrival of 
photons. 

The task can be formulated in very simple terms: If z is the detector position, at what 
time t is q a (t) = z a ? In other words: How to invert the equation (49), namely z a = 
q a + c (p a /p) t, to get tl Several comments are in order here. First, we are promoting t to 
the category of a q-number, while demoting q(t) to a given, external, parameter. This is 
the very task to accomplish to define a time operator. Second, the evolution equations that 
we have to invert to obtain the time of arrival is the set (49) of three equations, one per 
component, depending on a unique parameter t. To be compatible, they have to satisfy the 
constraint: 

l> a (z) = t abc (q b - z b ) p c = (50) 

To quantize this constrained system we borrow from the method of Dirac [50]. Classically, 
the constraint guarantees that the orbital angular momentum of the particle is z A p, so 
that z is a point of its trajectory. The total Hamiltonian formed by adding the constraints 
to the original one is: 

H T (z) = c^ + VaK(z) (51) 

where q a and p a are the dynamical variables to become operators after quantization, the 
/i's are Lagrange multipliers, and z is an external vector parameter corresponding to the 
detection position. The system is first class: 

{L a (z), L b (z)} = e abc L c (z), {L a (z), H T (z)} = e abc fi b L c (z) (52) 

where {, } indicate Poisson brackets as this is still classical dynamics. The evolution of the 
constraints does not produce additional (secondary) constraints. Hence, the Hamiltonian 
(51) is enough to account for the evolution of the constrained system. 

When quantizing the system, the vanishing of the constraint translates into a kind of 
subsidiary condition: 

(!„(*)$). (p)=0 (53) 
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The set of vectors of the Hilbert space that satisfy the above equation form the 

subspace H z of the states that could, eventually, be detected at z. On the other hand, using 
(27), Eq. (53) can be written as: 

I d 

e ahc ^]2\p\Y,ti(pA) {ij; z b } p c (pA|$) = (54) 

whose solution is 

(pA|$,) = e-** $(|p|,A,z) (55) 

Note that the functions $ may depend on the modulus of the momentum (but not on 
its direction), on the helicity and, possibly, on the detection point z. However, this last 
dependence has to be switched off to maintain translational invariance. 

We now proceed to invert Eq. (49). First of all we write down the action of the position 
operator on the detected subspace: 

q a (pA|$,)= e"** { z * + i^-JL} $(] P |,A) (56) 

Hence, with z a a parameter and i(z) an operator, Eq. (49) reads: 

z a = q a + c fj t(z), e~ ipz % fj $(|p|, A) + c fr t(z) e~ tpz $(|p|, A) = (57) 

\p\ IpI u\p\ \p\ 

This equation has to be valid whatever the function <3> chosen. This serves to define £(): It 
is precisely the operator that transforms (57) into an identity in H z , that is: 

i(z) w -i e- ipz -f- e ipz , t(z) = -i e- ipz ^- -f- \p\ e ipz (58) 
o\p\ \p\ o\p\ 

The symbol ~ at the left indicates equal up to operator ordering, something that we fix by 
the condition that t be maximally symmetric in the integration by parts with the measure 
d 3 p of the \p\) basis. This produces the operator defined at the right hand side, that we 
shall use as the time of arrival operator in what follows. It depends parametrically on z in 
as much the same way as the operators depend parametrically on t. Incidentally, this makes 
us recall that we are working in the Heisenberg picture at t — 0. It is straightforward to 
show that, when some time t has elapsed, the time operator shifts to t(z,t ) = t(z) — t , 
and that the arrival occurs at a time t z such that t(z,t z ) = 0. 

The eigenfunctions (p\\tz) of i(z) obtained by solving the eigenvalue equation 
t(z) (p\\tz) = t(p\\tz) are proportional to expi(Ht — pz)/\p\. Due to the fact that 
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lim|p|_>o(|p| (pMtty) 7^ 0' being |p| = the lower bound of the Hamiltonian, the op- 
erator (58) can not be self-adjoint. This is similar to the case of the radial momentum 
p r = 2r _1 (<9/<9r)r in three space dimensions (but notice that this does not prevent us from 
using H = pi /2m + 1(1 + l)/2mr 2 as Hamiltonian; what is done in this case is to restrict 
the behavior of the wave function at the boundary r = 0). The time operator commutes 
with the helicity so, taking into account that (p\\q\') = 5\\i(27r)~ 3 / 2 exp (—ipo), we could 
write 

(p X'\t X;z) = (p X'\ \-e lHt \z X) \t A; z) = \-e lHt \z A) (59) 

H H 

This notation is highly symbolic mainly due to the fact that z is just an external parameter 
belonging to the experimental set-up, the observer's will, etc, so that there is nothing like 
z among the properties of the particle. However, the above expression is correct if one 
considers that \z A) is the eigenket of the particle's position operator with eigenvalue z. 
Writing the detected states (55) in the same form may through some light on the meaning 
of the notation: 

(p\\Q x ) = (p\\Q(H, A))|z A) (60) 

This is the effect of the subsidiary condition (53): it projects on the detector position z, 
keeping only that part of the state that is in s-wave relative to z (hence the $ dependence on 
\p\ alone). The lack of completeness this produces on the time operator, and the associated 
interpretation, was discussed with some detail in ref. [47] for the relativistic massive spinless 
particle. I shall repeat here the two main results tailored to the massless, helicity ±1, photon 
case: 

1. The time eigenstates are not orthogonal: 

Y,[d 3 p (tX; z\pX'){pX'\t' A"; z) 

A' J 

= ^wr^^'^wv^kn (61) 

2. The basis is complete only within H z . In other words, the projection over states orthog- 
onal to the detector position is excluded from the decomposition of the identity in terms of 
time eigenstates: 

/+°° 2ir 
dt (p X\ t X'; z) ( t A'; z\ p" A") = 5\ \» -- 2 6(\p\ - |p»|) (p X\z X ) (z X |p" A" ) 
-°° \P\ 

/+oo 
dt (p A| t A'; z) ( t A'; z\ $ x ) = (p A| $ x ) V <t> z e H z (62) 

A' -°° 
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We recall that due to the non-orthogonality of the states, the spectral decomposition of 
the operator 



defines a positive operator valued measure only (not a projector valued one) that can be used 
to give the probability that the arrival of the state ip at the detector occurs in the time t as 
P^(t; z) = J2x |( t A; z\ ip )| 2 . This really is a probability density, something not reflected 
on the notation for the sake of simplicity. Note also that, due to the lack of completeness, 
the probability of eventually arriving at z (in any time) P^{z) — J dt P^(t;z) < 1, and 
may be zero for states ip ^ Ti z . Finally, the mean value of the time of arrival at z of a 
particle in the state ip is t^(z) = ( ip \i(z)\ ip )/P^(z) ) this excludes counterfactuals (the 
case Pjf,(z) = 0) as it should. 



Photons are eventually detected through their interaction with matter. Gauge invari- 
ance singles out the minimal coupling of the potentials to the currents J^(x,t) of additively 
conserved quantum numbers (the electric charge in this case). The coupling to Pauli like 
currents (e. g. d u ipa^ip) would never produce the finite amplitudes for absorbing or emit- 
ting soft photons observed experimentally. These facts and their implications have been 
thoroughly analyzed by S. Weinberg [31]and other authors. Of course, by means of a canon- 
ical transformation [51, 52] the minimal coupling interaction can be cast into a multipolar 
form. We will take into account this structure of the interactions when discussing the detec- 
tion of photons. Assume for simplicity the case of broad-band photodetectors for which the 
counting rates are given by the energy density of the field (ip\E^(x, t) ■ E{x, t) \ip). Then, the 
probability that the state ip be localized (in the plain sense of being detected by a detector) 
in a neighborhood A around (x,t) is given in terms of the fraction of the total energy that 
is within A. 

Absolute localization of quantum mechanical systems [5] -that is, the condition that the 
probability of finding the system out of some finite volume vanish- is such a strong condition 
that it violates causality [53, 54]. In other words, any free particle initially confined in a finite 
volume, continues in it forever, or immediately spreads to infinity. This result applies to free 
relativistic and non relativistic particles, to complex systems, in the presence of interactions, 




A 



(63) 



Conclusions 
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etc. Surprisingly the only requirement for this is that the system Hamiltonian be bounded 
from below. This prompts the question of what is the maximal degree of localization to be 
expected for a particle. An important clue [55] is that - even if the localization outside the 
finite volume is not absolute, but exponentially bounded tails are permitted - the probability 
spreads out to infinity faster than with any finite propagation speed. 

The limits to photon localization have been strengthened recently [26] based on simple 
physical requirements to be satisfied by one photon states. Basically, these are: a) That 
(pX\if)) can be given a probabilistic interpretation (something that we discussed below (16)) 
and b) That the Hamiltonian be bounded from below. Then, the Paley- Wiener Theorem 
VIII [56] says that the fall-off of the photon wavefunction (qX\ip) as \q\ — > oo is slower 
than exp(— a|q| r ), where a, r are positive constants and r < 1. As the physical requirements 
noted above apply to all types of particles, the same occurs to the limit of almost exponential 
localization: it applies to all kind of particles. This puts photons at the same level than 
the other particles in what refers to localization. A recent analysis of spontaneous emission 
from excited atoms [27] has shown the possibility of producing entangled atom-photon states 
where the photon wave packets have Gaussian tails. This explicit breaking of the barrier of 
exponential localization is a product of the entangled final state. 

These results shall likely find their application in the field of quantum information, and 
shall promote new developments in quantum optics. Some necessary tools like good po- 
sition [1] and time of arrival operators and their associated probabilistic interpretations 
are provided in this paper. We are completing a detailed analysis of the application of 
these tools to the tunneling through photonic band gaps, to HOM [57] interferometry 
and entanglement [27, 58], and to the superluminal propagation detected in several ex- 
periments [59, 60, 61]. 
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